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Abstract. A microscopic framework of nuclear energy density functionals is re- 
' viewed, which establishes a direct relation between low-energy QCD and nuclear 

structure, synthesizing effective field theory methods and principles of density 
functional theory. Guided by two closely related features of QCD in the low- 
energy limit: a) in-medium changes of vacuum condensates, and b) spontaneous 
breaking of chiral symmetry; a relativistic energy density functional is developed 
. and applied in studies of ground-state properties of spherical and deformed nuclei. 

33 
S3 . 

1 Introduction 

' Modern nuclear structure is rapidly evolving from studies of nuclei close to the /^-stability line 

towards regions of exotic short-lived nuclei and systems at the nucleon drip-lines. A wealth of 
new experimental results has also prompted important qualitative and quantitative advances in 
nuclear structure theory. The principal objective is to build a consistent microscopic framework 
that will provide a unified description of bulk properties, nuclear excitations and reactions. 

Among the microscopic approaches to the nuclear many-body problem, the framework of 
nuclear energy density functionals (NEDF) is the only one that can be applied over the whole 
chart of nuclides, from relatively light systems to superheavy nuclei, and from the valley of 
/^-stability to the particle drip-lines. NEDFs are not necessarily related to any realistic nucleon- 
^ . nucleon interaction, but rather represent global functionals of nucleon densities and currents. 

With a small set of universal parameters adjusted to data, EDFs have achieved a high level 
of accuracy in the description of structure properties. Nevertheless, one expects that new data 



00 
O 



OO 

O 

> 

• rH 

X 

on exotic nuclei with extreme isospin values will present interesting challenges for the NEDF 
framework. 

Even though completely phenomenological EDFs have been successfully employed in exten- 
sive structure studies, it would clearly be desirable to have a fully microscopic foundation for 
a universal EDF framework. Eventually, the description of the nuclear many-body problem, 
including both extended nuclear matter and finite nuclei, must be linked to and constrained 
by the underlying theory of strong interactions. In these lectures we will review a novel pro- 
gram [112131 1| , which integrates effective field theory (EFT) methods and principles of density 
functional theory (DFT), to bridge the gap between features of low-energy non-perturbative 
quantum chromodynamics (QCD) and nuclear phenomenology Guided by two important el- 
ements that establish connections with chiral dynamics and the symmetry breaking pattern 
of low-energy QCD: a) strong scalar and vector nucleon fields related to in-medium changes 
of QCD vacuum condensates; b) the long- and intermediate-range interactions generated by 
one-and two-pion exchange, as derived from in-medium chiral perturbation theory; a relativis- 
tic nuclear EDF has been developed which includes second-order gradient corrections to the 
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local density approximation. Translated into a relativistic point-coupling model with density- 
dependent vertices, this framework has been applied in calculations of ground-state properties 
over a broad range of spherical and deformed nuclei. 

Sec. [2] introduces elementary concepts of density functional theory, and in Sec. [3] we discuss 
some basic issues in the construction of nuclear EDFs. Elements of low-energy QCD are reviewed 
in Sec. |4j In-medium chiral EFT is applied in the calculation of the nuclear matter equation of 
state in Sec.[5l In Sec.[6]we develop a relativistic EDF based on chiral pion-nucleon dynamics, 
and present some illustrative calculations of nuclear ground-state properties. A brief summary 
is included in Sec. [7] 



2 Basic Concepts of Density Functional Theory 

Density Functional Theory (DFT) is one of the most popular and successful "ab initio" ap- 
proaches to the structure of quantum many-body systems (atoms, molecules, solids). Probably 
no other method achieves comparable accuracy at the same computational cost. The basic 
concept is that the ground-state properties of a stationary many-body system can be repre- 
sented in terms of the ground-state density alone. Since the density n(r) is a function of only 
three spatial coordinates, rather than the 3N coordinates of the N-body wave function, DFT is 
computationally feasible even for large systems. 

In this section we introduce elementary concepts of DFT, that will also be useful in con- 
sidering nuclear energy density functionals. Detailed expositions are available in numerous 
review articles, textbooks and collection of research papers. We refer the reader, for instance, 
to [5I6I7I8I9] , The basis of DFT are the Hohenberg-Kohn theorem and the Kohn-Sham scheme, 
which will be introduced in the next two sections. We will also briefly discuss the approximations 
to the central ingredient of the DFT approach: the exchange-correlation energy functional. 

2.1 The Hohenberg-Kohn Theorem 

In ground-state DFT one considers a system of N interacting particles described by the Hamil- 
tonian: 

N 9 N N N 

i— 1 i—l i— 1 j^i 

Let \P be the N-body wave function and n(r) the corresponding particle density: 

n(r) — N J d 3 r 2 J d 3 r 3 ...J d 3 r N !P*(r, r 2 , . . . , r N )&(r, r 2 , . . . , r N ) . (2) 

The Hohenberg-Kohn (HK) theorem [TIT] can be summarized in the following statements: 

— The nondegenerate ground-state (GS) wave function is a unique functional of the GS density: 

^ (ri,r 2 ,...,rjv) = #[no(r)] . (3) 
As a consequence, the GS expectation value of any observable O is also a functional of no(r): 

O = O[n ] = {*[n ]\d\¥[n ]). (4) 

- The GS energy: 

E = E[n ] = {¥[no]\H\9[no]) , (5) 

and the GS density no (r) of a system characterized by the external potential vo (r) can be 
obtained from a variational principle which involves only the density, i.e. the inequality: 

E =E Vo [n ] <E Vo [n] , (6) 

holds for all other densities n(r). 
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— There exists a functional F[n] such that the energy functional can be written as: 

E V0 [n] = F[n] + J d 3 r v (r)n(r) . (7) 

The functional F[n] is universal in the sense that, for a given particle-particle interaction 
(for instance the Coulomb interaction), it does not depend on the potential vq(t), i.e. it has 
the same functional form for all systems. 

The proof is based on the Rayleigh-Ritz variational principle |5|7|10j . The formal definition of 
the Hohenberg-Kohn functional F[n] then reads: 

F[n] = T[n] + W[n] = (W[n]\T\$[n]) + (&[n]\W\W[n]) , (8) 

where \P[n] is the N-body wave function which minimizes the expectation value of T + W. The 
GS density uniquely determines the external potential of the system and thus, as the kinetic 
energy and particle-particle interaction are universal functionals, the entire Hamiltonian and 
all physical properties of the system. However, although it establishes the variational character 
of the energy functional, the HK theorem gives no practical guide to the construction of the 
universal functional, and the explicit density dependence of F[n] is not known. In order to make 
practical use of DFT, one needs reliable approximations for T[n] and W[n], 



2.2 The Kohn-Sham Scheme 

Most practical applications of DFT use the effective single-particle Kohn-Sham (KS) equations, 
which are introduced for an auxiliary system of N non-interacting particles described by the 
Hamiltonian: 

H s =f + V s . (9) 
The HK theorem states that there exists a unique energy functional 

E s [n] = T s [n] + [ d 3 r »»(r)n(r) , (10) 



for which the variational equation yields the exact ground-state density n s (r) corresponding to 
H s , and T s [n] is the universal kinetic energy functional of the non-interacting system. The KS 
scheme is based on the following assertion: for any interacting system, there exists a unique 
local single-particle potential v s (r), such that the exact GS density of the interacting system 
equals the GS density of the auxiliary non-interacting system: 

N 



n(r)=n s (r)=^|^(r)| 2 , (11) 



expressed in terms of the lowest N single-particle orbitals - solutions of the Kohn-Sham equa- 
tions: 

[-V 2 /2m + v s (r)] &(r) = e^v) . (12) 

The uniqueness of v s (r) follows from the HK theorem, and the single-particle orbitals are unique 
functionals of the density: <j>i(r) — (/>j ( [n] ; r) . The HK functional is partitioned in the following 
way: 

F[n] = T s [n] + U[n] + E xc [n] , (13) 
where U[n] is the Hartree functional: 

U[n] =~ / d 3 r ( dVn(r)w(r,r>(r') , (14) 
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and the last term is the exchange-correlation energy functional -E^n] which, by the definition 
Eq. (fT3|) . includes everything else - all the many-body effects. It can formally be defined by: 

E xc [n] = T[n] + W[n] - U[n] - T s [n] . (15) 

The corresponding effective Kohn-Sham potential therefore reads: 

v s [n](r) = v(r) + / d 3 r'w(r,r')n(r') + v xc [n](r) , (16) 



where the local exchange-correlation potential is defined by: 



SE xc [n] 

v xc [n r = . 17 

on(r) 



Since the effective potential depends on the density, the system of equations (|11[) . p^|) . and (JTC 
has to be solved self-consistently. This is the Kohn-Sham scheme of density functional theory 

By construction, therefore, the exact u s [n ](r) is the unique external potential which, for 
the non- interacting system, leads to the same physical GS density no(r) as that of N interact- 
ing particles in the external potential v(r). In the limit when one neglects E xc (v xc ), the KS 
equations reduce to the self-consistent Hartree equations. However, since it includes correlation 
effects, the full KS scheme goes beyond the Hartree-Fock approximation and, in addition, it 
has the advantage of being a local scheme. With an exact functional E xc , all many-body effects 
are in principle included. Thus the usefulness of the Kohn-Sham scheme crucially depends on 
our ability to construct accurate approximations to the exact exchange-correlation energy. 



2.3 Approximations for the Exchange-Correlation Energy 

The true exchange-correlation energy functional is universal, i.e. given the inter-particle inter- 
action, it has the same functional form for all systems. The construction of an expression for 
the unknown functional E xc is therefore the keystone of the Kohn-Sham DFT. One possible 
approach is to develop E xc from first principles by incorporating known exact constraints. An- 
other is empirical, a parametric ansatz is optimized by adjusting it to a set of data. Modern 
approximations for E xc typically combine both strategies. 

The simplest, but at the same time still useful, is the well known local density approximation 
(LDA). In LDA, the exchange-correlation energy functional reads: 

E xc A \ n \ = J o?rn{v)e^{n{v)), (18) 

where e x " 1 ' (n(r)) is the exchange-correlation energy per particle of a uniform infinite system 
of density n (electron gas, nuclear matter), and which can be obtained from an "ab initio" 
calculation. By its construction, the LDA is a priori expected to be valid for spatially slowly 
varying densities. Nevertheless, LDA proved to be a surprisingly accurate approximation also 
for systems with rapid density variations. This is a consequence of exact properties which the 
LDA inherits from the homogeneous system. 

Any realistic system is spatially inhomogeneous, and it is clearly useful to incorporate the 
information on the rate of density variation in E xc . In the generalized-gradient approximation 
(GGA) the exchange-correlation functional is written as: 

Eg GA [n] = J d 3 r /(n(r),Vn(r)) . (19) 

GGA is more accurate then LDA but, unlike the input e"™^ in LDA, the function / is not 
unique and, depending on the method of constructing /(n(r),Vre(r)), very different GGAs 
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can be obtained. For instance, while typical GGAs used in quantum chemistry [12] are based 
on parametric forms fitted to sets of selected molecules, exact constraints are incorporated in 
first-principles GGAs |13] used in most physics applications. 

The next rung on the ladder of approximations for E xc are "meta-GGA" functionals which, 
in addition to the density and its gradient, depend on the kinetic energy density of the occupied 
Kohn-Sham orbitals 



Unlike LDA and GGA, which are explicit functionals of the density alone, the meta-GGA 
functional 



also depends explicitly on the Kohn-Sham orbitals. Meta-GGA is the highest level of approx- 
imation which docs not include full non-locality. Approximations of even higher level of accu- 
racy incorporate increasingly complex ingredients, and include fully non-local functionals of the 
Kohn-Sham orbitals, occupied as well as unoccupied. 



3 Nuclear energy density functionals 

The most complete and accurate description of structure phenomena in heavy nuclei is currently 
provided by self-consistent mean-field (SCMF) models. A variety of ground-state properties and 
collective excitations, not only in medium-heavy and heavy stable nuclei, but also in regions of 
exotic nuclei far from the line of /3-stability and close to the nucleoli drip- lines, have been suc- 
cessfully described with mean-field models based on the Gogny interaction, the Skyrme energy 
functional, and the relativistic meson-exchange effective Lagrangian [14115] . In the mean-field 
approximation the dynamics of the nuclear many-body system is represented by independent 
nucleons moving in a self-consistent potential, which corresponds to the actual density distri- 
bution of a given nucleus. Important advantages of using the mean- field framework include the 
use of global effective nuclear interactions, the treatment of arbitrarily heavy systems including 
superheavy nuclei, and the intuitive picture of intrinsic shapes. 

The SCMF approach to nuclear structure is analogous to Kohn-Sham density functional 
theory. As we have shown in the previous section, DFT enables a description of a quantum 
many-body problem in terms of a universal energy density functional, and nuclear mean-field 
models approximate the exact energy functional, which includes all higher-order correlations, 
with powers and gradients of ground-state nucleon densities and currents. Although it models 
the effective in-medium interaction between nucleons, a general density functional is not nec- 
essarily related to any given NN potential. By employing global parameter sets, adjusted to 
reproduce empirical properties of symmetric and asymmetric nuclear matter, and bulk prop- 
erties of simple, spherical and stable nuclei, the current generation of SCMF models has been 
applied in numerous studies of structure phenomena over the whole nuclear chart. 

The earliest applications of DFT in nuclear structure |16|17|18] used the zero-range density- 
dependent effective Skyrme interaction. The corresponding Skyrme functional can be written 
as the most general energy-density functional in isoscalar and isovector density, spin density, 
current, spin-current tensor, kinetic density, and kinetic spin density, respectively [14j : 




(20) 




(21) 





s o(r) 





j T (r) = i(V'-V)pr(r,r')| r: 
r T (r) = V-VV(r,r')l 



Mr) 
T T (r) 



|(V - V)<8>s T (r,r') 
V-V's T (r,r')l , 



(22) 
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where a denotes the spin, and r the isospin of the nucleon. Isoscalar (T = 0) quantities are sums 
(po = p n + Pp), whereas isovector (T = 1) densities correspond to proton-neutron differences 
(pi = p n — pp). The Skyrme functional contains systematically all possible bilinear terms in the 
local densities and currents of Eq. (|22|) up to second order in the derivatives, which are invariant 
with respect to parity, time-reversal, rotational, translational and isospin transformations: 

^ Sk = { C t Pt + c t P Pt^Pt + C T T pttt + C T JS + C^ J p T V ■ J T 

T=0,1 

C S T s 2 T + C$ s s T ■ As T + Cf s T ■ T T + C^ s ( V • s T ) 2 + C T f T + C^ j s T • V x j T } (23) 

where the coefficients Ct in the isoscalar and isovector channels can be either constants or 
explicitly depend on the nucleon density. 

Even though the functional can be derived from the ground-state expectation value of 
a zero-range momentum-dependent force introduced by Skyrme, in modern applications the 
functional is parameterized directly by fitting the coefficients to nuclear ground-state data, 
without reference to any NN interaction. Over the last thirty years more than hundred different 
Skyrme parameterizations have been adjusted and analyzed, and about ten of these parameter 
sets are still used in extensive studies of nuclear structure phenomena. This means, however, that 
often it is difficult to compare results obtained with different models, also because they include 
different subsets of terms from the general functional Eq. ([23)) . Ideally, model dependences could 
be removed by including all terms allowed by symmetries [19] . However, it seems that available 
data can only constrain a subset of parameters, and one needs additional criteria for selecting 
the optimal energy density functional form. 

One of the major goals of modern theoretical nuclear physics is, therefore, to build a universal 
nuclear energy density functional (NEDF) [20) . Universal in the sense that the same functional 
form is used for all nuclei, with the same set of parameters. This framework should then provide 
a reliable microscopic description of infinite nuclear and neutron matter, ground-state properties 
of all bound nuclei, low-energy vibrations, rotational spectra, small-amplitude vibrations, and 
large-amplitude adiabatic properties. In order to formulate a microscopic NEDF, however, 
one must be able to go beyond the mean-field approximation and systematically calculate 
the exchange-correlation energy functional E xc [p\. Many-body correlations must be included 
starting from the relevant active degrees of freedom at low-energies characteristic for nuclear 
binding. 

Nuclei are aggregates of nucleons (and mesons), which are in turn clusters of quarks and 
gluons. The decription of nuclear many-body systems must therefore ultimately be linked to 
and constrained by the underlying theory of strong interactions - quantum chromodynamics 
(QCD). The success of the non-relativistic and covariant SCMF approach to nuclear many-body 
dynamics |14|15j , and the recent application of chiral effective field theory to nucleon- nucleon 
scattering and the few-body problem [21] , point to a unified microscopic framework based 
on density functional theory (DFT) and effective field theories (EFT). As we will show in 
the following sections, within this framework a NEDF can be developed in a systematic way 
by using an EFT of low-energy in-medium nucleon-nucleon interactions to construct accurate 
approximations to the exact exchange-correlation energy functional. 

Effective field theories are low-energy approximations to more fundamental theories. Instead 
of trying to solve explicitly the underlying theory, low-energy processes are described in terms 
of active degrees of freedom. The low-energy behavior of the underlying theory must be known, 
and it must be built into the EFT. High-momentum states that characterize the unknown 
short-distance dynamics must be excluded, and the EFT is used to calculate observables in 
terms of an expansion in p/A (p denotes momenta or masses smaller than a certain momentum 
scale A). The resulting effective Lagrangian must include all terms that are compatible with 
the symmetries of the underlying theory and this means, in principle, an infinite number of 
terms. A working EFT must, therefore, be completed by specifying a method that allows to 
decide which terms of the effective Lagrangian contribute in a calculation up to a desired level 
of accuracy. 
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In the next section we will outline basic concepts of chiral EFT, a low-energy approximation 
to QCD. At low energies characteristic for nuclei, the active degrees of freedom of QCD are 
pions and nucleons, and their dynamics is governed by the chiral SU(2) x SU(2) symmetry and 
its spontaneous breakdown. 



4 Elements of Low-Energy QCD 

Quantum chromodynamics is the fundamental quantum field theory of strong interactions. 
At very short distance scales, r < 0.1 fm, corresponding to momentum transfers exceeding 
several GeV/c which probe space-time intervals deep inside the nucleon, QCD is a theory of 
weakly interacting quarks and gluons. At momentum scales considerably smaller than 1 GeV 
(corresponding to length scales r > 1 fm, typical for the physics of atomic nuclei), QCD is 
characterized by color confinement and a non-trivial vacuum: the ground state of QCD hosts 
strong condensates of quark-antiquark pairs and gluons. The associated spontaneous breaking of 
chiral symmetry implies the existence of pseudoscalar Goldstone bosons. For two quark flavours 
(u and d quarks), the Goldstone bosons are identified with the isospin triplet of pions. At low 
energy, pions interact weakly with one another and with any massive hadron. At the energy 
scale of nuclear binding, QCD is thus realized as an effective field theory of pions coupled to 
nucleons. 

The quark condensate (qq), i.e. the ground state expectation value of the scalar quark 
density, plays a particularly important role as an order parameter of spontaneously broken 
chiral symmetry Hadrons and nuclei represent excitations built on this condensed ground state. 
Changes of the condensate structure of the QCD vacuum in the presence of baryonic matter 
are a source of strong fields felt by the nucleons which, together with the pion-exchange forces, 
govern nuclear dynamics. 

In this section we briefly outline some basic concepts of QCD in the low-energy two-flavour 
(Nf = 2) sector. For a more detailed introduction we refer the reader to Refs. [22123124] . 

4.1 The QCD Lagrangian 

In the low-energy sector of QCD the lightest u and d quarks form a flavour Nf = 2 (isospin) 
doublet with "bare" quark masses of less than 10 MeV. The flavour (and colour) components 
of the quarks are collected in the Dirac fields tp(x) = (u(x), d(x)) T . The QCD Lagrangian 

Cqcd = Tp(ij„D» - m)V - ^Tr(G^Gn, (24) 
includes the gluonic field strength tensor 

= d„Ap a - d v Ap a + gf abc A<^A b J c , (25) 

the covariant derivative 

D^ = d^- %gAp a , (26) 

and the generators t a = A a /2 (a = 1,...,8) of SU(3) co i OI group. The 2x2 matrix m = 
diag(m u , rrid) contains the light quark masses. Since quarks are not observed as asymptotically 
free states, the quark mass depends on the momentum scale at which it is probed. At a scale 
of the order of 1 GeV: m u — 4 ± 2 MeV, rrid — 8 ± 4 MeV [35] . The strange quark is more than 
an order of magnitude heavier (m s ~ 150 MeV), and it will not be considered here. 

The interaction between quarks and gluons is independent of the quark flavour. The coupling 
strength a s = g 2 /An depends on the momentum scale Q. At high energies, a s is small and QCD 
is characterized by asymptotic freedom. In this regime it is possible to apply perturbation 
theory with an expansion in a s . For Q < 1 GeV, a s increases sharply and the relevant degrees 
of freedom are no longer the elementary colored quarks and gluons. They are replaced by color 
neutral bound states: the hadrons. Color confinement characterizes QCD at low energies. 
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4.2 Chiral Symmetry 

Let us consider QCD in the limit of massless quarks, by setting m — in Eq. (|24|) . In this limit, 
the QCD Lagrangian has a global symmetry related to the conserved right- or left-handedness 
(chirality) of zero mass spin 1/2 particles. Introducing right- and left-handed quark fields, 

lfo,i = i(l±7feM (27) 

we observe that separate global unitary transformations 

fa - exp[i0%?£] fa, fa - exp[i9 a L ^] fa, (28) 

with r a (a = 1, 2, 3) the generators of (isospin) SU(2), leave Cqcd invariant in the limit m — > 

£qcd fah^fa + fai^D^fa (29) 

This is the chiral SU(2) r x SU (2) ^ symmetry of QCD. It implies six conserved Noether currents, 
J R,a = $Rl t *T a ,/2fa> and j£ a = ■ipLj fi T a /2ip L , with 3^3% = d^j£ = 0. It is also useful to 
introduce the vector and axial currents, 

V£ = JR.a + J£,a = H*ty > = J ^ - JL,a = Witty ■ (30) 

Their corresponding charges, 

Ql = J d 3 x ^{x)ty{x) , QA = J d 3 x ^{ X ) lb ty{ x ) , (31) 

are, likewise, generators of SU{2) x SU{2). 

Chiral symmetry is explicitly broken by the finite quark masses, because the mass term in 
the Lagrangian mixes left- and right-handed fields: 

C m = -ipmip = -(ipRmfa + Tp L miJjR) . (32) 

As a consequence, the charge operators Eq. (I31[) are, in general, no longer time independent. 

4.3 Spontaneous Chiral Symmetry Breaking 

A (continuous) symmetry is said to be spontaneously broken or hidden, if the ground state 
of the system is not invariant under the full symmetry of the Lagrangian. In the case of the 
QCD Lagrangian with m = 0, the chiral SU (2) x SU(2) symmetry is spontaneously broken: the 
ground state (vacuum) of QCD has lost part of the symmetry of the Lagrangian. It is symmetric 
only under the subgroup SU(2)y generated by the vector charges Q v (isospin symmetry). 

Evidence for the spontaneous breaking of chiral symmetry is found in the spectrum of 
physical hadrons. If the QCD ground state were symmetric under chiral SU(2) x SU(2), it 
would be annihilated both by the vector and axial charge operators: Q^|0) = Q„ |0) = 0. This 
would then imply the existence of parity doublets in the hadron spectrum which, however, are 
not observed in nature. For instance, the vector (J" = 1~) p meson mass (m p ~ 0.77 GeV) is 
well separated from that of the axial vector (J 11 = 1 + ) a\ meson (m Q1 ~ 1.23 GeV). Likewise, 
the light pseudoscalar (J^ = 0~) mesons have masses much lower than the lightest scalar 
(J n = + ) mesons. 

Goldstone's theorem states that in a quantum field theory every spontaneously broken con- 
tinuous global symmetry leads to massless particles (Goldstone bosons) with the same quantum 
numbers as the generators of the broken symmetry. If G is the symmetry group of the La- 
grangian with rig generators, and H the subgroup with uh generators which leaves the ground 
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state invariant after spontaneous symmetry breaking, the total number of Goldstone bosons 
equals ng — n#. 

In the case of chiral symmetry, to each axial generator of SU(2)a, which does not 
annihilate the ground state, corresponds a massless Goldstone boson field with spin and 
negative parity. Namely, if Q„ |0) 7^ 0, there must be a physical state generated by the axial 
charge, \<p a ) = QaW)-> which is energetically degenerate with the vacuum. Let Hq be the QCD 
Hamiltonian (with massless quarks) which commutes with the axial charge. Setting the ground 
state energy equal to zero for convenience, we have H \(f) a ) — Q^H \Q) = 0. Evidently \<j> a ) 
represents three massless pseudoscalar bosons (for Nf = 2). They are identified with the pions: 

7T~, 7T°, and 7T + . 

Goldstone's theorem consequently implies non- vanishing matrix elements of the axial current 
between the vacuum and the pion (0|A M |7r), which must be proportional to the pion momentum: 

(OK^Mg)) = -iUeS^e-** , (33) 

and the constant of proportionality is the pion decay constant f v = 92.4 MeV, determined from 
the decay n + — ► + /i + f M 7. The divergence of Eq. ((33)) reads: 

(0\dMx)\TT b (q)) = -UtfSabe-*** = -Uml5 ab e-^ x . (34) 

The small pion mass (~ 140 MeV), as compared to hadronic scales, is directly linked to the 
partial conservation of the axial current (PC AC) Eq. ([34]) . The axial transformation is thus an 
approximate symmetry of QCD, and becomes exact in the limit m — > 0. 



4.4 The Chiral Condensate 

Spontaneous symmetry breaking is related to a scalar operator whose nonvanishing vacuum 
expectation value plays the role of the order parameter of the symmetry breaking. The QCD 
ground state contains scalar quark-antiquark pairs, and the corresponding expectation value 
(0 1 -0-0 1 0} is called the chiral (or quark) condensate: 

(tpip) = (uu) + (dd) = -Tr lim (0\Ti/>{x)tp(y)\0) , (35) 

where T denotes the time-ordered product. The relation between spontaneous chiral symmetry 
breaking and the non-vanishing chiral condensate can be derived from the (equal-time) com- 
mutator between the pseudoscalar operator P a (x) = r tjj(x)^5T a 'ip(x) and the axial charge of 
eq. (f3Tj) : [Q^ , P b ] = —5 a b'4>ip- Taking the ground state expectation value of this commutator, we 
notice that Q„ |0) 7^ is consistent with (ipifi) 7^ 0. 

The Goldstone boson (pion) acquires its physical mass through the explicit breaking of 
chiral symmetry by the finite quark masses m Uj( j. The pion mass is related to the u- and 
d- quark masses by the Gell-Mann, Oakes, Renner (GOR) relation [26 : 

ml = --p{m u + m d )(qq) + 0{ml d ). (36) 

Isospin symmetry (Q„ |0) = 0) has been used in the definition: (qq) = (uu) ~ (dd). Neglecting 
terms of order d , identifying / with the empirical pion decay constant f v = 92.4 MeV, and 
inserting m u + m d ~ 12 MeV 25 27], one obtains 

(qq) ~ -(240 MeV) 3 ~ -1.8 /to" 3 . (37) 

This condensate presents a measure of the degree of spontaneous chiral symmetry breaking. 
The non-zero pion mass, on the other hand, reflects the explicit symmetry breaking by the 
small quark masses, with ~ m q . 
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Fig. 1. Chiral condensate (in units of its vacuum value) as a function of temperature and baryon 
density (po — 0.16 fm -3 is the saturation density of nuclear matter). 

The strength of the scalar condensate Eq. (|37[) is more than an order of magnitude larger 
than the nuclear matter density at saturation: po — 0.16 fm -3 . Fig. [T] displays the tempera- 
ture and baryon-density dependence of the chiral condensate in units of its vacuum value. At 
relatively low temperatures which are characteristic for low-energy nuclear physics, we notice 
an approximately linear dependence on baryon density up to p ~ po, and even beyond. The 
strength of the quark condensate at normal nuclear densities is reduced by about one third from 
its vacuum value. The temperature dependence is much weaker, at least up to T < 100 MeV, 
far above the temperatures that can be sustained by ordinary nuclei. The density-dependent 
changes of the condensate structure in the presence of baryonic matter are a source of strong 
scalar and vector fields experienced by the nucleons. At nuclear matter saturation density sev- 
eral hundred MeV of scalar attraction are compensated by an almost equal amount of vector 
repulsion. The sum of the condensate nucleon fields almost vanishes in infinite homogeneous 
nuclear matter. 



4.5 Chiral Effective Field Theory 

Chiral EFT is a low-energy approximation to QCD. In the hadronic, low-energy phase, the 
active degrees of freedom of QCD are not elementary quarks and gluons, but rather mesons 
and baryons. The EFT describes the active light particles as collective degrees of freedom, and 
the heavy particles are treated as static sources. The QCD Lagrangian is replaced by an effective 
Lagrangian which includes all relevant symmetries of the underlying fundamental theory. Chiral 
EFT can then be used to calculate physical processes in terms of an expansion in p/A, where 
p denotes momenta or masses smaller than a certain momentum scale A. 

In the meson sector (baryon number B = 0) the elementary quarks and gluons are replaced 
by Goldstone bosons. In the case of interest here, for two quark flavours (Nf — 2), the pseu- 
doscalar Goldstone pions are represented by a unitary 2x2 matrix field U(x) S SU(2) which 
collects the three isospin components ir a (x) (a = 1, 2, 3), and transforms under chiral rotations 
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ipR —> RijjR and %jj L — » Lip L , as: 

U^LUR^ with U(x) = exp[iT a <p a (x)} , (38) 

and <fi a = 7T //, where / is the pion decay constant in the chiral limit. Small corrections to 
/, of the order of the u- and d- quark masses, yield the physical pion decay constant f n = 
f + 0(m u , d ) = 92.4 MeV. 

Pions interact weakly at low energy: if \tt) = Q A \0) is a massless state with H\n) = 0, then a 
state \n n ) = (Q A ) n \0) with n pions is also massless since the axial charges Q A all commute with 
the full Hamiltonian H. Interactions between pions must therefore vanish at zero momentum 
and in the chiral limit. 

It is easy to show that 

Tr^Ud^] -» ^[Ld^UR^Rd^tf] = Tr^Ud^} (39) 

is invariant under chiral transformations. The QCD Lagrangian of Eq. (|24jl is replaced by an 
effective Lagrangian which involves an expansion in the field U(x) and its derivatives: 

Cqcd ^ C eSf (U,dU,d 2 U,...). (40) 

Momenta and derivatives are equivalent for this expansion. Pions do not interact unless they 
carry non-zero four-momentum, so the low-energy expansion of Eq. (|40[) is an ordering in 
powers of d^U. Only even numbers of derivatives are allowed, because the Lagrangian has to 
be a Lorentz sacalar: 

£ e// = 4 0) + 4 2) + 4 4) + - (4i) 

However, since U is unitary, can only be a constant. The expansion converges if the 
momenta are much smaller than a characteristic scale. In the case of chiral perturbation theory, 
this scale is A x = Airf^ — 1161 MeV, which provides a natural separation between "light" and 
"heavy" degrees of freedom. The leading term, which is also know as the "non-linear sigma 
model" , involves two derivatives: 

£ (2) = i-Tr^^d^U}. (42) 

A completely symmetric Lagrangian corresponds to massless Goldstonc bosons, whereas the 
physical pions are massive. This is taken into account by adding a symmetry breaking term 
linear in the quark mass matrix m: 

4 2 ) = ^Tr[d^d»U] + ^-B Q Tr[m(U + [/*)] . (43) 

The symmetry breaking mass term is small and it can be treated perturbatively, together with 
the power series in momentum. The index 2 denotes either two derivatives or one quark mass 
term, and this Lagrangian contains already all possible terms up to second order. At fourth 
order, the terms permitted by symmetries are : 

4 4) - j(Tr[d^d^U}) 2 + l -^Tr[d^d,U}Tr[d^d u U] , (44) 

and the constants h,h (following the standard notation of ref. [28] ) must be determined by 
experiment. The fourth order Lagrangian contains also symmetry breaking terms, with 
additional constants U not determined by chiral symmetry. 

An expansion of to terms quadratic in the pion field yields the pion Lagrangian plus a 
constant vacuum contribution: 

4 2 ) = (m u + m d )f 2 B Q + ^d^ a d^ a - X -(m u + m d )B Q n 2 a + 0(tt 4 ) . (45) 
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The first term in Eq. (|45j) corresponds to the constant shift of the vacuum energy density by 
the non-zero quark masses, and the pion mass term is identified as m 2 . = (m u + rrid)Bo. Using 
the GOR relation Eq. (J36j> , one finds -f 2 B = (qq). 

In addition to the expansion of the effective Lagrangian in terms with an increasing number 
of derivatives and quark mass terms, Chiral Perturbation Theory (ChPT) is defined by a sys- 
tematic method of estimating the importance of diagrams generated by the interaction terms - 
the power counting scheme. Diagrams are classified according to the power of the generic vari- 
able Q (three- momentum or energy of the pion, or the pion mass m ff ), and the small expansion 
parameter is Q/A x . 

In order to apply ChPT to the nucleon sector, pion-nucleon interactions must consistently 
be included in the computational framework based on the low-energy expansion of the effective 
Lagrangian. The free-nucleon Lagrangian reads: 

=¥ N (yy»& t -M )¥ Ii (46) 

where the Dirac spinor \&n(x) = (p, n) T denotes the isospin-1/2 doublet of proton and neutron, 
and Mq is the nucleon mass in the chiral limit. The low-energy effective Lagrangian for pions 
interacting with a nucleon is obtained by replacing the pure meson Lagrangian by C e ff — 
+ + ... + C^ff , which also includes the nucleon field. The additional pion-nucleon term 
is expanded in powers of derivatives (momenta) and quark masses: 

C?„ + (47) 

and the structure of the irN effective Lagrangian at each order is again determined by chiral 
symmetry. The leading term, C^, reads: 

4w = <M^(d" - *V) + 9A1^ - M ]* N . (48) 

It includes the chiral covariant derivative with the vector coupling between the pions and the 
nucleon, and the axial-vector coupling. At this order the Lagrangian contains two parameters 
not determined by chiral symmetry: the nucleon mass Mq and the axial-vector coupling constant 
QA. The vector and axial vector fields can be expressed in terms of the chiral matrix £ = \f\J : 

^ = ^d^ + ^^) = -^e abc T a ir b d»Tr c + ... , (49) 

i(£t^_^t ) = __l 



^ = -(£t0M£_£0M£t ) = _ Tad ^ a + ... , (50) 



and expanded in powers of the pion field and its derivatives 

(2) 

At next-to-leading order, in the symmetry breaking quark mass term appears, and 
shifts the nucleon mass from its value in the chiral limit to the physical mass: 

M N =M +a N . (51) 



The "sigma term" cr/v represents the contribution from explicit chiral symmetry breaking to 
the nucleon mass 

<jn = rn„— ■ — — = (N\m u uu + rridddlN) , (52) 

q—u,d 

through the non-vanishing u- and c?-quark masses. Its empirical value cat = (45 ± 8) MeV, has 
been deduced from low-energy pion-nucleon scattering data [29j . 

The 7riV effective Lagrangian, expanded to second order in the pion field, reads: 

C% f = &n(V - M N )V N - ljj-* NW T& N ■ d»iT (53) 



J TT J TT 
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The microscopic origin of the empirical axial- vector coupling constant qa = 1.267 ± 0.003, 
determined from neutron beta decay (n — > pep), is in the substructure of the nucleon, that 
cannot be resolved at the level of a low-energy effective theory. 

For a detailed review of baryon chiral perturbation theory we refer the reader to Ref . [30j . In 
applications of in-medium ChPT to the nuclear many-body problem (cf. Sec.[S]), the calculations 
will be based on the leading terms in the effective Lagrangian, linear in the derivative 9 M 7r of 
the pion field. 



4.6 QCD Sum Rules 

The QCD sum rule approach [31] has been designed to interpolate between the perturbative 
(short-distance) and non-perturbative (large distances) energy sectors. It provides qualitative 
and quantitative information on hadron parameters (masses and coupling constants) from vac- 
uum expectation values of QCD operators. 

The basic quantity is a momentum-space correlation function (correlator) of color-singlet 
currents: 

n{q 2 ) = i f d 4 x e l "- x (017^7(20^(0)10) (54) 



where Jix) is a current composite operator of quark fields, T denotes the time-ordered product, 
and |0) is the physical non-perturbative vacuum. The singularities of the correlator correspond 
to physical excitations with the quantum numbers of the current J . The sum rule analysis 
proceeds by introducing the spectral function of the correlator defined by the dispersion relation 

n(q 2 )= [°°ds 2 P( ; } . , (55) 
Jo s - 1 - ie 

where 

O(qo)p(q 2 ) = -ImTTGz 2 ) = (2^) 3 V 5\p k - q) | (0|J(0) \k) | 2 (56) 

fc 

describes the spectrum of physical intermediate states. For a fixed three-momentum q, the 
spectral function measures the intensity at which energy is absorbed from the current at different 
frequencies. 

The second step in the sum rule approach is a direct calculation of the correlator using the 
operator product expansion (OPE), which provides a QCD approximation to the correlator 
that is applicable at large spacelike q 2 . The basic concept is the expansion of a time-ordered 
product of two local operators at short distances in terms of a complete set of local operators: 

TA(x)B(0) x =° J2 C n B (x)On(0). (57) 

n 

The c-number coefficients C^ B of the expansion are called Wilson coefficients. Therefore, in 
moment urn- space : 

n(q 2 ) = i I 'd 4 x e^ x {Q\TJ{x)j\Q)\Q) = £ C„(<z 2 )(0|O„|0) (58) 
J n 

The operators O n are built from quark and gluon fields. The operators of lowest dimension are 
m q qq and G" G£". At the next order one finds four-quark operators and combinations of quark 
and gluon fields. The right-hand side of Eq. (|58[) is an expansion in vacuum expectation values 
(condensates) of the operators O n . The OPE is in inverse powers of Q 2 = —q 2 , and converges 
rapidly for large spacelike values of Q 2 > 0. 

The OPE separates short-distance and long-distance physics. In QCD this separation is 
between perturbative physics (the coefficients) and non-perturbative physics (the condensates) . 
The Wilson coefficients decrease with inverse powers of Q 2 for each higher dimension operator, 
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and can be calculated in QCD perturbation theory. The non-perturbative physics is contained 
in the condensates. 

The fundamental relation of the QCD sum rule method is obtained by equating the dis- 
persion relation for the correlator 77 (q 2 ) and its OPE representation, each evaluated at large 
spacelike Q 2 : 



n( q 2 ) = - [ 

Jo 



, ImiT(s) 

as -= + polynomial 



Q 

= C (-Q 2 ) + d(-Q 2 )(m g OT) + C 2 (~Q 2 )(G; u G^) + ■■■ (59) 

In applications of QCD sum rules the convergence of both sides of Eq. (|59|) is improved by 
applying the Borel transform - a differential operation defined, for a given function f(Q 2 ), by 
the relation: 

with Q 2 /n = M 2 . The transform Bf(Q 2 ) = f{M 2 ) depends on the parameter M. - the Borel 
mass. Any simple polynomial in Q 2 is eliminated by the Borel transform, because it cannot 
survive the infinite number of differentiations, and an inverse power of Q 2 is replaced by an 
exponentially decreasing function of the Borel mass: 

= -T72^ S/M2 ■ ( 61 ) 
s + Q 2 M 2 



Therefore, the Borel transform applied to Eq. f|59[) simultaneously eliminates the subtraction 
terms accompanying the dispersion relation and any divergent polynomial from the OPE. 
Higher-order terms in the OPE, which contain inverse powers of Q 2 , are suppressed factorially 
by the Borel transform: 

1 _ 1 
(Q^T ~ {n-l)\(M 2 Y 



The correlator in spectral form can be evaluated by introducing a phcnomcnological model 
for the spectral function Eq. (|56|) . Maching the Borel transforms of the OPE and phenomeno- 
logical description of the correlator, yields the sum rules for each invariant function. The leading 
order result for the nuclcon mass in vacuum is the Ioffe's formula: 



-V.s =-^a<9?) + -" (63) 



For M. « 1 GeV, one finds Mjv ~ 1 GeV. Ioffe's formula is not very accurate and the contribu- 
tion from higher-dimensional condensates should also be included. Nevertheless, it demonstrates 
that the scale of the nucleon mass is largely determined by the scalar quark condensate, i.e. by 
the spontaneous chiral symmetry breaking mechanism. 

In section [S~4l we will show how QCD sum rules at finite densities can be used to relate the 
leading changes of the scalar quark condensate (qq) and of the quark density (q^q) at finite 
baryon density, with the scalar and vector self-energies of a nucleon in the nuclear medium. 



5 Chiral Dynamics and Nuclear Matter 

In this section we begin to approach one of the central problems of modern theoretical nuclear 
physics: how to establish a relationship between low-energy, non-perturbative QCD and the 
rich nuclear phenomenology, which includes both nuclear matter and finite nuclei. We will first 
focus on infinite nuclear matter and analyze the equations of state of symmetric and asymmetric 
matter, as well as pure neutron matter, in the framework of chiral effective field theory. It will 
be shown that nuclear binding and saturation arise primarily from chiral (pionic) fluctuations 
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Fig. 2. NN amplitude in chiral effective field theory: one-pion exchange, two-pion exchange (including 
A isobar intermediate states), and contact terms representing short-distance dynamics. 

in combination with Pauli blocking effects and three- nucleon (3N) interactions, superimposed 
on the condensate background fields and calculated according to the rules of in-medium chiral 
perturbation theory (ChPT). 

The key clement of chiral EFT is a separation of long- and short- distance dynamics and 
an ordering scheme in powers of small momenta. In a medium, the relevant quantity for this 
expansion is the Fermi momentum kf, which defines the nucleon density via the relation: 

At nuclear matter saturation density (po ~ 0.16 fm -3 ): kfo ~ 2m n , so the Fermi momentum 
and the pion mass are of comparable magnitude at the densities of interest. This means that 
pions should be included as explicit degrees of freedom in the description of nuclear many-body 
dynamics. The relevant observable is the energy per particle E(kf), i.e. the ratio of energy 
density and particle density, with the free nucleon mass M subtracted. 

Both kf and m w are small compared to the characteristic chiral scale, 47r/„. ~ 1.2 GeV. Con- 
sequently, the equation of state (EOS) of nuclear matter as given by ChPT will be represented 
as an expansion in powers of the Fermi momentum. The expansion coefficients are non-trivial 
functions of kf /m^, the dimensionless ratio of the two relevant scales inherent to the problem. 



5.1 In-medium Chiral Perturbation Theory 

In the chiral EFT framework nuclear dynamics is described by a low-energy effective represen- 
tation of QCD. The relevant degrees of freedom are pions and nucleons, and observables are 
calculated in chiral perturbation theory, as described in Sec. 14.51 The chiral effective Lagrangian 
Eq. (|48[) . generates the basic pion- nucleon coupling terms. In the nuclear medium, however, all 
nucleon states below some Fermi momentum kf are occupied. Instead of the "empty" vacuum 
|0), the ground state is the filled Fermi sea \4>o), and the nucleon propagator in coordinate space 
changes to: 

S$(x-y) = (Q\T1>(x)$(y)\0) -> S F (x-y) = {<h\T^(x)$(y)\<h) . (65) 

In momentum space there is a very useful additive decomposition of the in-medium nucleon 
propagator. For a relativistic nucleon with four-momentum p^ = (po,p) it reads 

U + M ){ p 2_M2 + i£ - 2 ^P 2 - M 2 )e( Po )9(kf - |p |)j . (66) 

The second term is the medium insertion which accounts for the fact that the ground state of 
the system has changed from an "empty" vacuum to a filled Fermi sea of nucleons. Diagrams 
can then be organized systematically in the number of medium insertions, and an expansion is 
performed in leading inverse powers of the nucleon mass, consistent with the fc/-expansion. 

The NN amplitude is represented by the diagrams shown in Fig. [2] One- and two-pion 
exchange processes are treated explicitly. They govern the long-range interactions at distance 
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Fig. 3. One-pion exchange Fock diagram, and iterated one-pion exchange Hartree and Fock diagrams. 

o 

Fig. 4. Irreducible two-pion exchange diagrams. 

scales d> 1/kf fm, whereas short-range dynamics is not resolved at nuclear Fermi momentum 
scales, and can be subsumed in contact interaction terms. 

We notice that the two-pion exchange diagram in Fig. [5] includes not only nucleon interme- 
diate states, but also single and double virtual Z\(1232)-isobar excitations. In fact, in addition 
to kf and m w , a third relevant "small" scale is the mass difference 5M = Ma — Mjv ~ 0.3 GeV 
between the Z\(1232) and the nucleon. The strong spin-isospin transition from the nucleon to 
the A isobar must therefore be included as an additional ingredient in the in-medium ChPT 
calculation, so that the expansion coefficients of the nuclear matter equation of state become 
functions of both kf/m^ and m w /SM. 

Initial applications of in-medium ChPT to nuclear dynamics }32[33j have shown that basic 
properties of the equation of state for isospin-symmetric nuclear matter can already be repro- 
duced by 7riV-dynamics alone, without introducing virtual Z\-excitations. In the next section, 
therefore, we will first discuss nuclear matter properties calculated in a "one-parameter" chi- 
ral approach, which considers only 7riV-dynamics. In Sec. 15.31 we will show that the explicit 
inclusion of Z\(1232) degrees of freedom leads to a systematic improvement of the chiral EFT 
description of nuclear matter, especially in the isovector channel. 



5.2 Chiral "one-parameter" Nuclear Matter Equation of State 

The in-medium ChPT three-loop calculation of the energy per particle of symmetric nuclear 
matter has been carried out in Ref. [33], including contributions up to order 0(kf). These 
contributions are: the kinetic energy (contributing at order 0{kf)), chiral one-pion exchange 
(0(k 3 f )), iterated one-pion exchange (0(kf)), and irreducible two-pion exchange (O(fc^)). The 
corresponding closed-loop diagrams are shown in Fig. [3] (the one-pion exchange Fock diagram, 
and iterated one-pion exchange Hartree and Fock diagrams), and in Fig. |4] (irreducible two- 
pion exchange diagrams). The 27r Fock diagram has both a reducible part contributing to the 
iterated l7r exchange, and an irreducible two-pion exchange part [34] . 

The calculation involves one single momentum-space cut-off A which encodes dynamics at 
short distances not resolved explicitly in the effective low-energy theory. This high-momentum 
scale A is the only free parameter. Note that instead of a momentum-space cut-off, one could 
equivalently use dimensional regularization, remove divergent loop integrals, and replace them 
by adjustable NN contact terms which then parameterize the unresolved short-distance physics. 

The saturation of nuclear matter can simply be illustrated in the exact chiral limit: m n = 0. 
The saturation mechanism can already be demonstrated by truncating the one- and two-pion 
exchange diagrams at order O(kj). The energy per particle, E(kf) = E(kf)/A, of isospin 
symmetric nuclear matter is given in powers of the Fermi momentum kf, and can simply be 
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parameterized: 



lb 2 k 3 b 4 

— On, f n, f nj f 

E ^^k~ a ik +P ik- (67) 



The first term is the kinetic energy, and the result for a in the chiral limit is: 



V 47T / M N V 47T / 



- 1 



(68) 



where Mjy is the nucleon mass, and g^N — 9aMn/ fn is the empirical irN coupling constant: 
g nN = 13.2. The strongly attractive leading /c^-term is accompanied by the weakly repulsive 
one-pion exchange Fock term. The /^-contribution to E{kf) would lead to collapse, were it not 
for the stabilizing fc^-term controlled by the coefficient: 



P = 10 O 4 (4 ^ + 237 - 24 In 2) - | ^ 13.55 . (69) 

The two-pion exchange between nucleons generates the attraction proportional to fc 3 in the 
energy per particle. The Pauli blocking of intermediate nucleon states in these 2n exchange 
processes produces the repulsive kj term. The parameter-free ChPT value (3 = 13.55 in the 
chiral limit, is very close to the value (~ 12.2) that one obtains by adjusting the simple equation 
of state Eq. (|6T[) to empirical saturation properties (binding energy and saturation density) [33j . 
For values of the short-distance scale A between 0.5 and 0.6 GeV the EOS has a stable minimum 
of E(kf) in the empirical range of density and binding energy. 

The full 3-loop chiral dynamics result for E(p) in symmetric nuclear matter with nucleon 
density p(kf) — 2fc 3 /(37r 2 ), using = 135 MeV, and including all calculated terms up to order 
O(kj), is shown in the upper left panel of Fig. [5l The single-parameter of the calculation - the 

cut-off scale A = 0.65 GeV, has been adjusted to the value Eq = —15.3 MeV at equilibrium, 
and the resulting energy per particle E(p) has a minimum with this value of Eq at the density 
po = 0.178 fm -3 (corresponding to a Fermi momentum of kfo = 272.7 MeV = 1.382 fm _1 ). 
The predicted compression modulus: 



d^Ejkf) 



255 MeV, (70) 



kf—kf 



is in quantitative agreement with the "empirical" value 250 ± 25 MeV 35 36J. It is also inter- 
esting to analyze the various contributions to the equilibrium value Eq = —15.26 MeV [33]. 
Its decomposition into contributions from the kinetic energy and the three classes of diagrams 
(l7r-exchange Fock, iterated l7r-exchange, and irreducible 27r-exchange) reads: Eq = (23.40 + 
18.24—68.35 + 11.45) MeV. The nuclear matter binding energy can also be decomposed into con- 
tributions ordered in chiral powers 0{k u f ), {v = 2,3,4,5): Eq = (23.75- 154.54+124.61-9.08) 
MeV. 

The simplest way to extend the in-medium ChPT calculation to isospin asymmetric matter 
is to use the same cut-off for all isospin channels. The following substitution: 

*(fc/-|p|) - ^0(*p-|pI) + ^0(*»-|pI). (71) 

in the in-medium nucleon propagator Eq. (|66p . takes into account the fact that the Fermi seas 
of protons and neutrons are no longer equally filled. k p and k n denote the Fermi momenta of 
protons and neutrons, respectively. Choosing k pn — fc/(l + S) 1 ^ 3 (with S a small parameter) 
the nucleon density p — p p + p n — (k 3 + k^J/Sw 2 = 2kj/3ir 2 remains constant. The expansion 
of the energy per particle of isospin asymmetric nuclear matter: 

E as (k p , k n ) = E{k f ) + S 2 A(k f ) + • • • , (72) 
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Fig. 5. Energy per particle of symmetric nuclear matter as function of the nucleon density, deter- 
mined by chiral one- and two-pion exchange up to order 0(k%) (upper left panel). The corresponding 
asymmetry energy as function of the nucleon density (upper right), and the density dependence of the 
energy per particle of pure neutron matter (lower left) . The momentum dependence of the real part of 
the single-nucleon potential at nuclear saturation Fermi momentum kfo = 272.7 MeV, is shown in the 
lower right panel. 



around the symmetry line (k p — k n , or (5 = 0) defines the asymmetry energy A(kf). Note that 
the parameter S is equal to (p„ — p p ) / (p n +p p ), or (N — Z)/ (N + Z) . In the upper right panel of 
Fig. [5] we display the density dependence of A(kf) as determined by 7riV-dynamics up to O(kj) 
and three-loop order. At the saturation point kfo = 272.7 MeV: Ao = A(kfo) = 33.8 MeV, in 
very good agreement with the empirical value: Aq = (30 ± 4) MeV. Extrapolation to higher 
density works roughly up to p ~ 1.5 po. At still higher densities one observes an unphysical 
downward bending of the asymmetry energy curve, indicating the limit of validity of the chiral 
expansion scheme restricted to pion-nucleon dynamics only, and with only a single momentum 
cut-off A = 0.65 GeV parameterizing the unresolved short-range dynamics. 

The extreme of asymmetric nuclear matter is pure neutron matter. This system is unbound 
and its energy per particle increases monotonically with neutron density. To calculate the energy 
per particle, E n (k n ) of neutron matter in the framework of in-medium ChPT, it is sufficient to 
make the substitution 

*(*/-|p|) - ^0(*n-|PD. (73) 

in the nucleon propagator Eq. (|66p . Here k n denotes the Fermi momentum of the neutrons, 
related to the neutron density by p n = k^/iii 2 . The resulting EOS of neutron matter is shown 
in the lower left panel of Fig. [5l The convex shape of the curve is generic and does not change 
much with the cut-off A. When compared with sophisticated many-body calculations of neutron 
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matter EOS, one finds a rough agreement up to neutron densities of about p n — 0.25 fm -3 . 
At higher densities the chiral ttN neutron equation of state becomes unrealistic because of the 
downward bending, inherited from the asymmetry energy A(kf). 

The in-medium three-loop calculation of the energy per particle defines the (momentum de- 
pendent) self-energy of a single nucleon in nuclear matter up to two-loop order. The momentum- 
and density-dependent (complex- valued) single-particle potential of nucleons in isospin symmet- 
ric nuclear matter has been calculated in Ref. [37] . The momentum dependence of the real part 
U(p,kfo) of the single-particle potential evaluated up to order O(fcl), at nuclear saturation 
Fermi momentum kfo — 272.7 MeV, is shown in the lower right panel of Fig. [5] For a nucleon 
at rest, one finds C/(0, fc/o) = —53.2 MeV, in agreement with the empirical optical-model po- 
tential Uo « —52 MeV, deduced by extrapolation from elastic nucleon-nucleon scattering data. 
The momentum dependence of U(p, fc/o) translates into an effective nucleon mass M^(p): 

! ' d {T kin (p) + U(p,k f0 )} . (74) 



m n(p) P d P 

The strong momentum dependence with the up- and down-bending shown in Fig. produces 
a negative slope of U(p, kfo) at p = kfo, resulting in a far too large nucleon mass at the Fermi 
surface. As we will show in the next section, a significant improvement of the momentum depen- 
dence of U(p, kfo) is obtained by explicitly including diagrams involving virtual Z\-excitations. 

The realistic values for several nuclear matter properties (binding energy and saturation 
density, compressibility, asymmetry energy, depth of the single- nucleon potential), have been 
obtained with only one adjustable scale parameter - the momentum cut-off A. These results 
demonstrate that the nuclear binding and saturation mechanism based on the combination 
of chiral one- and two-pion exchange, with short-distance dynamics encoded in a single high- 
momentum scale A, presents a very good starting point for a description of the nuclear many- 
body problem based on chiral dynamics, at least at low nucleon densities p < 0.2 fm~ 3 . 



5.3 The Role of Virtual ZX(1232) Excitations 

The chiral approach to nuclear matter has been extended and improved in Ref. [38) , by system- 
atically including contributions from two-pion exchange with single and double virtual Z\(1232)- 
isobar excitations. The spin-3/2 isospin-3/2 Z\(1232)-resonance is the most prominent feature 
of low-energy 7rV-scattering. Two-pion exchange between nucleons with excitation of virtual 
Z\-isobars generates most of the empirical strong isoscalar central NN-attraction [39], which 
in phcnomenological one-boson exchange models is often simulated by a fictitious scalar "a- 
meson" exchange. In addition, as we have already emphasized in Sec. 15.11 the delta-nucleon 
mass splitting SM — 293 MeV is of the same size as the Fermi momentum kfo — 2m^ at 
nuclear matter saturation density, and therefore pions and Z\-isobars should both be treated as 
explicit degrees of freedom. 

In Ref. [38] a calculation of isospin-symmetric and isospin- asymmetric nuclear matter EOS 
has been carried out, including all effects from 2-7r-exchange with virtual Z\-excitation up to 
three-loop order in the energy density. The leading contributions from 27r-exchange with vir- 
tual Z\-excitation to the energy per particle are generically of fifth power in the small momenta 
(kf, rriir, SM). The effects from irreducible 27r-exchange (cf. Sec. I5.2|) are of the same order. How- 
ever, since the 7r7VZ\-coupling constant is about twice as large as the 7rW-coupling constant 
{g-KNA = ig-nN I '\/2), the Z\-driven 27r-exchange effects should dominate. The contributions to 
the energy per particle can be classified as two-body terms and three-body terms, with the lat- 
ter interpreted as Pauli-blocking effects on the two-body terms imposed by the filled Fermi-sea 
of nucleons. 

Fig. [6] displays the relevant one-loop triangle, box and crossed box two-pion exchange di- 
agrams with single and double Z\(1232)-isobar excitations. By closing the two open nucleon 
lines to either two rings or one ring one gets (in diagrammatic representation) the Hartree or 
Fock contribution to the energy density of nuclear matter. The NN potential in momentum 
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Fig. 6. One-loop two-pion exchange diagrams with single and double Z\(1232)-isobar excitations. 



space involves pion-loop diagrams which are in general ultra-violet divergent and require regu- 
larization (and renormalization). This procedure results in a few subtraction constants which 
encode unresolved short-distance NN-dynamics. The associated & 3 and fc^-terms in the energy 
per particle 



m NN \k f )=B 3 



B f 



(75) 



are then adjusted to some empirical property of nuclear matter (e.g. the binding energy of ~ 16 
MeV). £?3 and B§ are chosen dimensionless, and Mn = 939 MeV is the nucleon mass. Eq. (175)) 
replaces the cut-off regularization used in Sec. 15.21 and, in addition, also takes into account the 
effects from p 2 -dependent contact terms. 

The resulting EOS for symmetric nuclear matter displays binding and saturation in a wide 
interval of values of the two adjustable parameters B% and B§. However, a strong repulsive 
p 2 -term from the three-body Hartree diagram leads to a too low saturation density po and a 
too high nuclear matter compressibility: K > 350 MeV. This problem has been solved in a 
minimal way by introducing an attractive three-body contact term, 



m NNN Hk f )=B 6 -fa . (76) 

1V1 N 

The three parameters: B 3 and B 5 of the two-body contact term Eq. {75]), and B 6 of the 
three-body term Eq. (|76[) . determine not only the isospin-symmetric energy per particle E{kf) 
at three-loop order, but also the the momentum dependence of the single-particle potential 
U(p, kf). With the following choice of the parameters: B3 = —7.99, B§ = 0, and Bq = —31.3 
[38] . the minimum of the saturation curve E(kf) is fixed to the value Eq = —16.0 MeV, the 
predicted value of the saturation density is po = 0.157 fm -3 , corresponding to a Fermi momen- 
tum of kfo = 261.6 MeV = 1.326 fm -1 , and the nuclear matter compressibility equals K = 304 
MeV, a somewhat high but still acceptable value. The density dependence of the energy per 
particle is shown in the upper left panel of Fig. [71 in comparison with the previous EOS cal- 
culated without the inclusion of virtual zi-excitations (cf. Sec. I5.2[) . In that case the saturation 
density po — 0.178 fm -3 was somewhat too high, but the compressibility K — 255 MeV had a 
better value. The more pronounced increase of the solid curve with the density p results from 
the inclusion of higher order terms in the ^/-expansion. 

The solid curve in the lower right panel in Fig. [7J represents the real part of the single- 
particle potential U(p,kfo) at saturation Fermi momentum kfQ = 261.6 MeV, as a function 
of the nucleon momentum p, whereas the dashed curve shows the single-particle potential 
calculated with the restriction to only nucleon intermediate states in 27r-exchange diagrams. 
We observe that with the chiral 7rVZ\-dynamics included, the real single-particle potential 
U(p,kfo) increases monotonically with the nucleon momentum p. The unphysical downward 
bending above p = 180 MeV displayed by the dashed curve, is no longer present. The slope 
at the Fermi surface p = kfo translates into an effective nucleon mass of M*(kfo) — 0.88M, a 
realistic value compared to M*(kfo) ~ 3M obtained in the previous calculation. The potential 
is now somewhat deeper: U(0, fc/o) = — 78.2MeV, compared to U(0, fc/o) = —53.2 MeV for the 
case without virtual Z\-excitations, but still realistic. 
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Fig. 7. Energy per particle of symmetric nuclear matter (upper left), and the asymmetry energy (upper 
right) as functions of the nucleon density. Energy per particle of pure neutron matter (lower left) as 
a function of neutron density, and the momentum dependence of the real part of the single-nucleon 
potential at nuclear saturation (lower right). The dashed curves refer to the results described in Sec. 15. 21 
with only pions and nucleons as active degrees of freedom. The solid curves include the contributions 
from two-pion exchange with single and double virtual Z\(1232)-isobar excitations. 



The equations of state of pure neutron matter are compared in the lower left panel of Fig. 
[7] With the inclusion of chiral 7riVZ\-dynamics, the short-range contribution reads: 



E^\k n ) = B n>3 ^+B n<5 -^ 



M 2 



(77) 



where B n s and i? n ,5 are two new subtraction constants. Note that the Pauli-exclusion principle 
forbids a three-neutron contact-interaction. The (dashed) solid curve represents the energy per 
particle E n (k n ) of pure neutron matter as a function of the neutron density p n = k^/Sir 2 , 
calculated (without) with the inclusion of chiral 7r./VZ\-dynamics. The short-range parameters 
B n ,3 — —0.95 and B n ^ — —3.58 have been adjusted to the value of the asymmetry energy at 
saturation density A(kfo) — 34 MeV (cf. upper right panel). The important result here is that 
the unphysical downward bending of E n (k n ) (displayed by the dashed curve) disappears after 
the inclusion of irNA diagrams. It has also been noted that up to p n » 0.16 fm -3 , the chiral 
ttNA result for E n (k n ) is similar to the standard EOS obtained with sophisticated many-body 
calculations based on free-space NN-interactions plus adjustable 3N-forces, whereas it becomes 
stiffcr at higher densities. Of course one should not expect the present model to work at Fermi 
momenta larger than k n w 350 MeV. 

The improved isospin channel of the EOS is also manifest in the density dependence of the 
asymmetry energy A(kf) (upper right panel of Fig. [7]), defined by Eq. (|72]l. The important 
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feature is that the inclusion of chiral 7rA^Z\-dynamics (solid curve) eliminates the unrealistic 
downward bending of the asymmetry A(kf) at higher densities p > 0.2 fm -3 (dashed curve). 
The value at saturation density po = 0.157 fm -3 : A(fc/o) = 34.0 MeV, is consistent with the 
empirical values of the asymmetry energy: (30 ± 4) MeV. 

In these last two sections we have shown that a microscopic approach to the nuclear many- 
body problem, based on chiral EFT, and which includes pions and Z\-isobars as explicit degrees 
of freedom in in-medium ChPT calculations, successfully describes the equations of state of 
isospin-symmetric and isospin-asymmetric nuclear matter. In fact, a nuclear equation of state 
with realistic saturation properties (binding energy, saturation density, compression modulus) 
can already be generated by the simplest version of the chiral framework, including only the 
contributions of one- and iterated one-pion exchange up to order 0{kj) and parameterizing all 
unresolved short-range dynamics with a simple momentum cut-off. In the calculation restricted 
to only nucleon intermediate states, i.e. without virtual Z\-excitations, the asymmetry energy 
at the saturation point is well reproduced and neutron matter is predicted to be unbound, but 
the curves for both the asymmetry energy and the neutron matter equation of state display an 
unphysical downward bending at higher densities. The real part of the single-nucleon potential 
has the correct depth at momentum p = 0, but the strong momentum dependence results in 
a far too large nucleon mass at the Fermi surface. These deficiencies of the in-medium ChPT 
calculation are cured by the explicit inclusion of contributions from two-pion exchange with 
single and double virtual Zi(1232)-isobar excitations. With two parameters per isospin channel 
for the unresolved short-range physics, the expansion is carried out up to three-loop order and 
up to fifth order in small scales (kj, m^, 8M). The chiral 7r7VZ\-dynamics considerably improves 
both the momentum dependence of the single particle potential (realistic value for the effective 
mass), and the isovector properties of nuclear matter (density dependence of the asymmetry 
energy and neutron matter EOS). 



5.4 In-medium QCD Sum Rules 

Initial applications of the microscopic approach based on chiral ttN dynamics (not yet including 
the Z\(1232)) to bulk and single-nucleon properties of finite nuclei |l|2j . have shown that chiral 
(two-pion exchange) fluctuations play a prominent role in nuclear binding, but do not produce 
the strong spin-orbit nuclear force responsible for the large energy spacings between spin-orbit 
partner single-nucleon states. The large effective spin-orbit potential in finite nuclei is generated 
by additional strong isoscalar scalar and vector nucleon self-energies, induced by changes of the 
quark condensate and the quark density in the presence of baryonic matter. The magnitude 
of these condensate background nucleon self-energies can be estimated using leading-order in- 
medium QCD sum rules. 

Applications of QCD sum rules at finite baryon density p proceed analogously to the vacuum 
case (cf. Sec. I4.6j) . except that now operators such as q^q also have non-vanishing ground state 
expectation values: (q^q) p = 3p/2. In addition, the vacuum condensates vary with increasing 
density. Thus in-medium QCD sum rules relate the changes of the scalar quark condensate and 
the quark density at finite baryon density, with the isoscalar scalar and vector self-energies of 
a nucleon in the nuclear medium. In leading order which should be valid at densities below 

and around saturated nuclear matter, the condensate part, £g*\ of the scalar self-energy is 
expressed in terms of the density dependent chiral condensate as follows [40141142] : 

^(0) 8tt 2 8tt 2 a N 

where p s = is the nucleon scalar density. The chiral vacuum condensate (qq)o as a measure 
of spontaneous chiral symmetry breaking in QCD is given in Eq. (|37[) . The difference between 
the vacuum condensate (qq)o and the one at finite density involves the nucleon sigma term, ctat ~ 
(N\m q qq\N) , to this order. The Borel mass scale M. w lGeV roughly separates perturbative 
and non-perturbative domains in the QCD sum rule analysis. 
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To the same order in the condensates with lowest dimension, the resulting time component 
of the isoscalar vector self-energy is 

(o) 64^ 2 f 32vr 2 

It reflects the repulsive density-density correlations associated with the time component of the 
quark baryon current, qj^q. Note that around nuclear matter saturation density where p a ~ p, 
the ratio 

.(0) 

(80) 



4(m„ + m d ) p 

is approximately equal to —1 for typical values of the nucleon sigma term <tn and the current 
quark masses m M d (as an example, take ctat ~ 50 MeV and m u + rrid — 12 MeV at a scale of 1 
GeV). 

Identifying the free nucleon mass at p — with the Ioffe's formula Eq. (|63|) . M — — xj2(f<7)o, 
one finds 

sf{p)=M* { p)-M = -^i Ps (81) 



and 



with the quark masses m Ul d taken at a renormalization scale p ~ M. ~ Anf^ ~ 1 GeV. Eq. 1(81 
implies (identifying M with the free nucleon mass): 



(0). . 4(to m + m d )M 

(P) = -775 P > (82) 



4 0) ~ -350 MeV — ^— ^ . (83) 
50 MeV po 



Evidently, the leading-order in-medium change of the chiral condensate is a source of a strong, 
attractive scalar field which acts on the nucleon in such a way as to reduce its mass in nuclear 
matter by more than 1/3 of its vacuum value. 

We note that the QCD sum rule estimates implied by Eqs. (JHTJ) , f|82jl and by the ratio 
Eq. (|80|) are not very accurate at a quantitative level. The leading-order Ioffe formula on which 
Eq. (f8Tj) relies, should be corrected by contributions from condensates of higher dimension 

which are not well under control. The estimated error in the ratio Eg /Ey ~ — 1 is about 
20%, given the uncertainties in the values of <tn and m u + ma- Nevertheless, the constraint 
implied by Eq. (|80p will be very useful as a starting point in constructing the relativistic energy 
density functional for finite nuclei. 



6 Relativistic Energy Density Functional Based on Chiral EFT 

Starting from the successful in-medium ChPT calculation of homogeneous nuclear matter ex- 
plored in the previous section, we will develop a relativistic energy density functional (EDF) 
and apply it to studies of the structure of finite nuclei. The construction of this EDF is based 
on the following conjectures |l|2|3j : 

1. The nuclear ground state is characterized by strong scalar and vector mean fields which 
have their origin in the in-medium changes of the scalar quark condensate and of the quark 
density. 

2. Nuclear binding and saturation arise primarily from chiral (pionic) fluctuations in combina- 
tion with Pauli blocking effects and three-nucleon (3N) interactions, superimposed on the 
condensate background fields. 
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6.1 The Nuclear Energy Density Functional 

From Sec. [2] where the basics of density functional theory (DFT) were introduced, we recall 
that the energy functional is commonly decomposed into three separate terms: 

F[p] = T km [p] + E H [p] + E xc [p] , (84) 

where Thin is the kinetic energy of the non-interacting N-particle system, Eh is a Hartree energy, 
and E xc denotes the exchange-correlation energy which, by definition, contains everything else. 
In the formulation of the relativistic EDF we assume that the large scalar and vector mean fields 
that have their origin in the in-medium changes of the chiral condensate and of the quark den- 
sity, determine the Hartree energy functional En[p], whereas the chiral (pionic) fluctuations 
including one- and two-pion exchange with single and double virtual Z\(1232)-isobar excita- 
tions plus Pauli blocking effects, are incorporated in the exchange-correlation energy functional 
E xc [p}. 

The density distribution and the energy of the nuclear ground state are obtained from 
self-consistent solutions of the relativistic generalizations of linear single-nucleon Kohn-Sham 
equations. In order to derive those equations it is useful to construct a model with point- 
coupling NN interactions and density-dependent vertices, designed in such way to reproduce 
the detailed density dependence of the nucleon self-energies resulting from Ejj[p] + E xc [p\. For 
a two component system of protons and neutrons we start from a relativistic Lagrangian which 
includes isoscalar-scalar (S), isoscalar- vector (V), isovector-scalar (TS) and isovector- vector 
(TV) effective four-fermion interaction vertices with density-dependent coupling strengths. The 
minimal Lagrangian density reads: 

•C = -Cfrcc + ^i n t + ^int ^coul > (85) 

with the four terms defined as follows: 

Arcc = i>{il^ - M N )i> , (86) 

-~ G TS {p){i> r t V)- \ Gtv(p)$ t 7m V) • $ ?M , (87) 

C cm = eA^^p-^ - - A F^F^ , (89) 

where ip is the Dirac field of the nucleon with its two isospin components (p and n). Vectors 
in isospin space are denoted by arrows. In addition to the free nucleon Lagrangian £f ree and 

the interaction terms contained in c\^ t , when applied to finite nuclei the model must include 
the coupling £ om of the protons to the electromagnetic field with — d^A v — d u A^, and 

(2) 

a derivative (surface) term £ in t- One could, of course, construct additional derivative terms 

(2) 

in £j nt , further generalized to include density-dependent strength parameters. However, there 



appears to be no need in practical applications to go beyond the simplest ansatz (|88l) with a 
constant D$- Note that we do not introduce explicit spin-orbit terms. They emerge naturally 
from the Lorentz scalar and vector self-energies generated by the relativistic density functional. 

The variational principle SC/Sip = applied to the Lagrangian Eq. (|85|) leads to the 
self-consistent single-nucleon Dirac equations, the relativistic analogue of the (non-relativistic) 
Kohn-Sham equations. The nuclear dynamics produced by chiral (pionic) fluctuations in the 
medium is now encoded in the density dependence of the interaction vertices. In the framework 
of relativistic density functional theory [7143144) , the density-dependent couplings are functions 
of the 4-current : 

f = = pu» , (90) 
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where is the 4- velocity defined as (1 — v 2 ) _1 / 2 (l, v). We work in the rest-frame of the nuclear 
system with v = 0. 

The couplings Gi{p) (i — S,V 7 TS,TV) are decomposed as follows: 

Gi(p) =G\ 0) +G<?\p) (bri = S,V) 
and Gi(p) = Gf\p) (for i = TS,TV) . (91) 

The density-independent parts arise from strong isoscalar scalar and vector background 

fields, whereas the density-dependent parts G^(p) are generated by one- and two-pion ex- 
change dynamics. It is assumed that only pionic processes contribute to the isovector channels. 

The relativistic density functional describing the ground-state energy of the system can be 
re-written as a sum of four distinct terms: 



E a [p] = E lrcc [p] + Eg [p] + E cou} [p] + E^ [p] , (92) 



with 



Efreelp] = J d 3 r (^[-ry ■ V + M N ]i!\<f> Q ) , (93) 
Eg[f>] = \J<? r {(<Aol4 O) (^) 2 |0o) + (MG^i^flfo)} , (94) 
EAP) = l J d 3 r {{^ol^H^^I^ + ^olG^^)^^) 2 !^) 

+(0o 1 41 (/3) (V 5 r V) 2 |0o) + (0o|4t(/3)(^ ^>)>o) 

- (4>o\DP[V(H)f\^)} , (95) 

1 f ,3 , ,u_ .7.1 + ^3. 



= 2j dr (0o|^— ^ 7 ,W#o> , (96) 

where |0o) denotes the nuclear ground state. Here £f ree is the energy of the free (relativistic) 
nucleons including their rest mass. Eg is a Hartree-type contribution representing strong scalar 
and vector background condensate fields, and E v is the part of the energy generated by chiral 
7r./VZ\-dynamics, including a derivative (surface) term. Minimization of the ground-state energy, 
represented in terms of a set of auxiliary Dirac spinors ipk , leads to the relativistic analogue of 
the Kohn-Sham equations. These single-nucleon Dirac equations are solved self-consistently in 
the "no-sea" approximation which omits the explicit contribution of negative-energy solutions 
of the relativistic equations to the densities and currents. This means that vacuum polarization 
effects are not taken into account explicitly, but rather included in the adjustable parameters 
of the theory [45|46j . 

The expressions for the isoscalar and isovector four-currents and scalar densities read: 

N 

h = EfcLi i>ki^k , j ^ = ^2 $kin t A , (97) 

k=l 
N 

PS =Y,k=l^k^k, Ps =^2$kTtp k , (98) 

fc=l 

where ifjf. are Dirac spinors and the sum runs over occupied positive-energy single-nucleon states. 
If we only consider systems with time-reversal symmetry in the ground-state, i.e. even-even 
nuclei, the space components of all currents vanish (j = 0) and, because of charge conservation, 
only the third component of isospin (rs = — 1 for neutrons and T3 = +1 for protons) contributes. 
The relevant combinations of densities are: 

p = (<h\&y°<l>\<h) = p p + p n , p3 = (0oI^t 37 V|0o) = p p - p n . (99) 

Ps = (<h\W>\<t>o) = Ps + Ps> PS3 = (<h\M\<l>o) = P P s-Ps- (100) 
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Minimization with respect to ?pk gives the single- nucleon Dirac equation: 

[-i-7 • V + M N + j S v + Jqt 3 S T v + lo^R + %s + T 3 S TS }^k = tfeV'fe , (101) 
with the self-energies: 

+ G^(p)}p + eA ^, (102) 

(P)P3, (103) 

+ G { ; ) (p)}p s + D { ; ) V 2 ps, (104) 

(p)psa, (105) 
1 { v\p) „ 2 , ggg(P) J 







Stv 


— "TV 




= [Gf 


Ets 


= G T J 




-if 



P + — ^7 — Ps- 



dp dp 

where A°(r) in Eq. (|102[> is the Coulomb potential. In addition to the usual contributions 
from the time components of the vector self-energies and the scalar potentials, we must also 
include the "rearrangement" term, Er, arising from the variation of the vertex functionals with 
respect to the nucleon fields in the density operator p. For a Lagrangian with density-dependent 
couplings, the inclusion of the rearrangement self-energies is essential in order to guarantee 

energy-momentum conservation and thcrmodynamical consistency p 2 -§^ f ~J = -| Yli=i T m (i-e. 
for the pressure equation derived from the thermodynamic definition and from the energy- 
momentum tensor). Using the single- nucleon Dirac equation and performing an integration by 
parts, the ground-state energy of a nucleus with A nucleons reads: 

£o = X> -\ I d 3 r{ [G<°> + gW(p)]pI + 41(p)P53 + [4 0) + G^(p)]p 2 + 

i i « 



fc=l 



Gtv(p) P3 + — q- — Ps P + — -Q p — PS3 P + 

dGP(p) p3+ dG^ plp+epchA0+D p psW2p ^ {m) 



dp dp 
where denotes the single- nucleon Kohn-Sham energies 



6.2 Low-energy QCD Constraints 



In Sec. 15.41 we have shown how the in-medium QCD sum rules relate the leading changes of the 
scalar quark condensate and of the quark density at finite baryon density, with the scalar and 
vector self-energies of a nucleon in the nuclear medium. Comparing Eqs. (|102|) and (|104p for 
the isoscalar vector and scalar potentials of the single- nucleon Dirac equations, with the Eqs. 
(|8Tj) and (|82f for the condensate background self-energies, respectively, the following estimates 
are obtained for the couplings of the nucleon to the background fields (the Hartree terms in 
the energy functional): 
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a N M. 
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,2 fl 



f 2 

J 7T 



-11 fin 



<7 N 



50 MeV 
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m d 



50 MeV 



(108) 
(109) 
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The many-body effects represented by the exchange-correlation density functional are ap- 
proximated by chiral 7T-/VZ\-dynamics, including Pauli blocking effects, fn the simplest DFT 
approach, the exchange-correlation energy for a finite system is determined in the local density 
approximation (LDA) from the exchange-correlation functional of the corresponding infinite 
homogeneous system, replacing the constant density p by the local density p(r) of the ac- 
tual inhomogeneous system. In our case the exchange-correlation terms of the nuclear density 
functional are determined within LDA by equating the corresponding self-energies in the single- 
nucleon Dirac equation (|I0I|) . with those arising from the in-medium chiral perturbation theory 
calculation of 7r7VZ\-dynamics in homogeneous isospin symmetric and asymmetric nuclear mat- 
ter (cf. Sec. GO): 

U(p= k f ,p) = ^ hPT (k f , P ) + S^{k f ,p) 
-U I (p=k f ,p)6 = £$f T (k f ,p) +S^ PT (k f ,p) , (110) 

where U(j>, kf) and Ui(p, kf) are the isoscalar and isovector momentum and density-dependent 
single-particle potentials, respectively, and S = (p n — p p )/(p n + p p ). In order to determine the 
density-dependent couplings of the exchange-correlation pieces generated by the point-coupling 
model, a polynomial fit up to order fcj is performed for the ChPT self-energies, and they are 
re-expressed in terms of the baryon density p = p p + p n and the isovector density p 3 — p p — p n : 

^S hPT ( fc /> P) = ( C sl + C s2P^ + C s3 pi + C siP ) p , (111) 

Z$ hPT (k f , p) = (c„i + c v2 p^ + c v3 pi + c v4 p) p , (112) 

ZrS PT (fc/, p) = (C ts l + C ts 2P^ + C ts3 p 3 + CtsAp) P3 , (113) 
£%y PT (kf,p) = (Ctvl + C tv 2P^ + C t ViP^ + CtviP) P3 , (114) 

Next these self-energies are identified with the corresponding contributions to the potentials 
Eqs. (|102[) - (|106[) in the point-coupling single-nucleon Dirac equation. The resulting expressions 
for the density-dependent couplings of the pionic fluctuation terms read ( small differences 
between p and p$ at nuclear matter densities are neglected here): 

G S (p) = C «l + C s2P^ + c s3 p^ + c s4 p , (115) 
G v ) (p) = 5 vi + c v2 p^ + c v3 pi + c v4 p , (116) 
G ts(p) = °tsi + c ts2 p^ + c ts3 p li + c tsi p , (117) 

Gtv(p) = Ctvl + Ctv2p^ + C tv3 p3 + CtV4P , (118) 



where the inclusion of the rearrangement term redefines the isoscalar-vector coefficients: 
Cvi = Coi, c V 2 = l/7(6c„ 2 - c s2 ), c v3 = l/4(3c„ 3 - c s3 ), and c v4 = l/3(2c„ 4 - c si ) 0. 

The coefficient of the derivative term in the equivalent point-coupling model (Eq. Q104[) ) 
can be determined from ChPT calculations for inhomogeneous nuclear matter. The isoscalar 
nuclear energy density emerging from chiral pion-nucleon dynamics has the form 38J: 

£(p,V P ) = P E(k f ) + (Vp) 2 F v (k f ) + ... . (119) 

-FV can be approximated by a constant in the relevant region of nuclear densities. The following 
relation between Ft? and the derivative term of the point-coupling model (see Eqs. (|88I95[) ) is 
then valid: 

- 2F V = D ( g ] . (120) 

The inclusion of derivative terms in the model Lagrangian and the determination of its strength 
parameters from ChPT calculations in inhomogeneous matter actually goes beyond the local 
density approximation. The term Eq. f88]). with the strength parameter given by Eq. ()120l) . 
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represents a second-order gradient correction to the LDA, i.e. the next-to-leading term in the 
gradient expansion of the exchange-correlation energy calculated by in-medium chiral pertur- 
bation theory. 

While bulk properties of infinite nuclear matter are useful for orientation, the large amount 
of nuclear observables provides a far more accurate data base that permits a fine tuning of the 
global parameters. The total number of adjustable parameters of the point-coupling model is 
seven, four of which are related to contact (counter) terms that appear in the ChPT treatment 
of nuclear matter. One parameter fixes a surface (derivative) term, and two more represent the 
strengths of scalar and vector Hartree fields. The values of the parameters are adjusted simul- 
taneously to properties of nuclear matter and to binding energies, charge radii and differences 
between neutron and proton radii of spherical nuclei, starting from the estimates for the cou- 
plings of the condensate background fields (Hartree term), and the constants in the expressions 
for the self-energies arising from chiral 7riVZ\-dynamics (exchange-correlation term). The result- 
ing optimal parameter set (FKVW) [3] is remarkably close to the anticipated QCD sum rule 
and ChPT values, with the exception of the two constants associated with three-body correla- 
tions, for which the fit to nuclear data systematically requires an attractive shift as compared 
to the ChPT calculation [38]. 

6.3 Nuclear Ground-State Properties 

The effective global FKVW interaction has been tested in self-consistent calculations of ground- 
state observables for spherical and deformed medium-heavy and heavy nuclei. The calcula- 
tions, including open-shell nuclei, are performed in the framework of the relativistic Hartree- 
Bogoliubov (RHB) model, a relativistic extension of the conventional Hartree-Fock-Bogoliubov 
method, that provides a basis for a consistent microscopic description of ground-state proper- 
ties, low-energy excited states, small- amplitude vibrations, and reliable extrapolations toward 
the drip lines |15j . The new microscopic FKVW interaction [3] has been employed in the 
particle-hole channel, in comparison with one of the most successful phenomenological effective 
meson-exchange relativistic mean-field interactions: DD-ME1 [47] • Pairing effects in nuclei are 
restricted to a narrow window of a few MeV around the Fermi level. Their scale is well separated 
from the scale of binding energies which are in the range of several hundred to thousand MeV, 
and thus pairing can be treated as a non-relativistic phenomenon. In most applications of the 
RHB model the pairing part of the well known and successful Gogny force [48 has been used 
in the particle-particle channel, and the same interaction is used in the illustrative examples 
included in this section. 

The isotopic dependence of the deviations (in percent) between the calculated binding ener- 
gies and the experimental values for even-A Pb nuclei, is plotted in the upper panel of Fig. [SI It 
is interesting to note that, although DD-ME1 and FKVW represent different physical models, 
they display a similar mass dependence of the calculated binding energies for the Pb isotopic 
chain. On a quantitative level the FKVW interaction produces better results, with the abso- 
lute deviations of the calculated masses below 0.1 % for A > 190. In lighter Pb isotopes one 
expects that the observed shape coexistence phenomena will have a pronounced effect on the 
measured masses. Because of the intrinsic isospin dependence of the effective single-nucleon 
spin-orbit potential, relativistic mean-field models naturally reproduce the anomalous charge 
isotope shifts. The well known example of the anomalous kink in the charge isotope shifts of 
Pb isotopes is illustrated in the lower panel of Fig. [8j The results of RHB calculations with the 
DD-ME1 and FKVW effective interactions are shown in comparison with experimental values. 
Both interactions reproduce in detail the A-dependence of the isotope shifts and the kink at 

208p| 3 

One of the principal advantages of using the relativistic framework lies in the fact that the 
effective single-nucleon spin-orbit potential arises naturally from the Dirac equation. The single- 
nucleon potential does not include any adjustable parameter for the spin-orbit interaction. In the 
FKVW model, in particular, the large effective spin-orbit potential in finite nuclei is generated 
by the strong scalar and vector condensate background fields of about equal magnitude and 
opposite sign, induced by changes of the QCD vacuum in the presence of baryonic matter. Fig. 
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Fig. 8. The deviations (in percent) of the calculated binding energies from the experimental values 
(upper panel) [49] . and the calculated charge isotope shifts in comparison with data [50], for the 
chain of even-A Pb isotopes. The charge isotope shifts are defined: Ar 2 h — r 2 h (A) — r 2 h ( 20S Pb) and 
Arl D = r 2 LD (A) - r 2 LD ( 208 Pb), where the liquid-drop estimate is r 2 LD (A) = '^rlA 2 ' 3 . 



12] displays the deviations (in percent) between the calculated and experimental values of the 
energy spacings between spin-orbit partner-states in a series of doubly closed-shell nuclei. The 
theoretical spin-orbit splittings have been calculated with the FKVW and DD-ME1 interactions. 
For the phenomenological DD-ME1 interaction the large scalar and vector nucleon self-energies 
which generate the spin-orbit potential, arise from the exchange of "sigma" and "omega" mesons 
with adjustable strength parameters. One notices that, even though the values calculated with 
DD-ME1 are already in very good agreement with experimental data, a further improvement 
is obtained with the FKVW interaction. This remarkable agreement indicates that the initial 
estimates for the condensate background couplings have been more realistic than anticipated, 
considering the uncertainties of lowest-order in-medium QCD sum rules. 

Deformed nuclei with N > Z present further important tests for nuclear structure models. 
Ground-state properties, in particular, are sensitive to the isovector channel of the effective 
interaction, to the spin-orbit term of the effective single-nucleon potentials and to the effective 
mass. The nuclear density functional constrained by low-energy QCD has been tested in the 
region 60 < Z < 80. Predictions of the RHB calculations for the total binding energies, charge 
radii and ground-states quadrupole deformations of even-Z isotopic chains have been compared 
with available data. With the FKVW effective interaction in the particle-hole channel, and 
pairing correlations described by the finite range Gogny D1S interaction, very good agreement 
with experimental values has been found not only for the binding energies and charge radii 
over the entire region of deformed nuclei, but excellent results have also been obtained for 
the ground-state quadrupole deformations. The level of agreement with data is illustrated in 
Fig. [TU] where, for the chains of Nd, Sm, Gd, Dy, Er, Yb, Hf, W, and Os isotopes, the calculated 
ground-state quadrupole deformation parameters /3a > proportional to the expectation value of 
the quadrupole operator (</>o|3z 2 — r 2 \<fio), are displayed in comparison with the empirical data 
extracted from B(E2) transitions. We notice that the RHB results reproduce not only the global 
trend of the data but also the saturation of quadrupole deformations for heavier isotopes. 
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Fig. 9. The deviations (in percent) between the theoretical and experimental values [51j . of the energy 
spacings between spin-orbit partner-states in doubly closed-shell nuclei. 



7 Concluding Remarks 

We have reviewed a novel microscopic framework of nuclear energy density functionals, which 
synthesizes effective field theory methods and principles of density functional theory, to establish 
a direct link between low-energy QCD and nuclear structure. At low-energies characteristic for 
nuclear binding, QCD is realized as a theory of pions coupled to nucleons. The basic concept of a 
low-energy EFT is the separation of scales: the long-range physics is treated explicitly (e.g. pion 
exchange) and short-distance interactions, that cannot be resolved at low-energy, are replaced 
by contact terms. The EFT building of a universal microscopic energy density functional allows 
error estimates to be made, and it also provides a power counting scheme which separates long- 
and short-distance dynamics. 

A relativistic nuclear energy density functional (EDF) has been introduced and constrained 
by two closely related features of QCD in the low-energy limit: a) in- medium changes of vac- 
uum condensates, and b) spontaneous chiral symmetry breaking. The changes of the chiral 
(quark) condensate and quark density in the presence of baryonic matter are sources of strong 
(attractive) scalar and (repulsive) vector fields experienced by nucleons in the nucleus. These 
fields produce Hartree mean-field nucleon potentials, and are at the origin of the large energy 
spacings between spin-orbit partner states in nuclei. 

The spontaneously broken chiral symmetry in QCD introduces pions as Goldstone bosons, 
and determines the pion-nucleon couplings. Starting from in-medium chiral perturbation theory 
calculation of homogeneous nuclear matter, the exchange-correlation part of the EDF is deduced 
from the long- and intermediate-range interactions generated by one- and two-pion exchange 
processes, with explicit inclusion of Zi(1232) degrees of freedom. Regularization dependent 
contributions to the energy density, calculated at three-loop level, are absorbed in contact 
interactions with constants representing unresolved short-distance dynamics. 

The EDF framework is realized as a relativistic point-coupling model with density-dependent 
interaction vertices. The solution of the system of self-consistent single-nucleon Dirac equations 
(relativistic analogues of Kohn-Sham equations) determines the nucleon densities which enter 
the energy functional. The construction of the density functional involves an expansion of nu- 
cleon self-energies in powers of the Fermi momentum up to and including terms of order 
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Fig. 10. Comparison between the RHB model (FKVW interaction plus Gogny pairing) predictions 
for the ground-state quadrupole deformation parameters of the Nd, Sm, Gd, Dy, Er, Yb, Hf, W, and 
Os isotopes, and experimental values [52] . 



or equivalently, 0(p 2 ) in the proton and neutron densities. The exchange-correlation energy 
functional, determined by chiral 7riVzi-dynamics in nuclear matter, is used in Kohn-Sham cal- 
culations of finite nuclei by employing a second-order gradient correction to the local density 
approximation. 

Excellent results have been obtained in studies of ground-state properties (binding energies, 
charge form factors, radii of proton and neutron distributions, deformations, spin-orbit split- 
tings) of spherical and deformed nuclei all over the chart of nuclides. The microscopic EDF 
reproduces data on the same level of comparison as the best fully phcnomcnological (non- 
relativistic and relativistic) self-consistent mean-field models. These results demonstrate that 
chiral EFT establishes a consistent microscopic framework in which both the isoscalar and 
isovector channels of a universal nuclear energy density functional can be formulated. Guided 
by the principles of density functional theory and based on chiral EFT, the construction of the 
nuclear EDF provides a quantitative ab initio description of the complex nuclear many-body 
problem starting from the fundamental theory of strong interactions. 
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